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BRIEF 


In  the  final  report  to  the  Office  of  Havel  Research,  on 
the  Nonr-669(l5)  pro gran  at  Syracuse  University,  is  enclosed 
one  complete  doctoral  dissertation  on  the  Thoneon  heat  in 
Metallic  System.  Since  this  work  comprised  a  major  effort 
for  the  past  three  years  at  Syracuse,  it  is  included  in  its 
entirety.  The  theoretical  study  end  experimental  data  from 
the  Thomson  heat  dissertation  is  to  he  presented  to  the  open 
literature  as  a  publication  vithin  the  following  year. 


ABSTRACT 


The  fundamental  concepts  of  the  Thomson  heat  and  the  basic  equa¬ 
tions  governing  the  heat  flovs  associated  with  its  generation  are  pre¬ 
sented.  From  considerations  of  these  equations ,  two  nev  techniques  for 
measuring  the  Thomson  coefficient  have  been  developed  and  examined  ex¬ 
perimentally.  The  tvo  nev  techniques  of  studying  a  sample  subjected  to 
thermal  and  electrical  gradients.  Included  a  non -equilibrium  approach, 
in  which  the  initial  rate  of  change  of  temperature  within  the  sa^le 
leads  to  a  value  of  the  Thomson  coefficient,  and  a  method  based  on  equi¬ 
librium  considerations.  The  latter  leads  to  a  very  simple,  straight¬ 
forward  equation  relating  the  Thomson  coefficient  to  the  change  in  equi¬ 
librium  temperature  associated  with  a  reversal  of  current  flow. 

An  experimental  apparatus  was  designed  and  utilised  to  examine  the 
validity  of  the  equations  which  were  developed.  The  tvo  techniques,  as 
well  as  a  standard  direct  measurement  technique,  which  involved  a  change 
in  current  to  compensate  for  the  temperature  change  resulting  from  a 
reversal  of  Thomson  heat,  have  been  applied  to  measurements  of  the 
Thomson  heat  of  nickel  and  lead.  The  measurements  were  made  in  the  tem¬ 
perature  range  of  UO°C,  where  very  little,  if  any,  reliable  data  is 
available.  The  results  obtained  with  these  tvo  nev  techniques  agree 
reasonably  well  with  the  data  obtained  utilising  the  standard  technique 
and  data  extrapolated  into  this  teaperature  range  from  previous  work. 
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I.  INTRODUCTION 


The  phenomena  of  the  zmoele  Ctrl  city  encompasses  those  processes  In 
which  there  is  generation  of  electrical  currents  by  means  of  imposed 
thermal  gradients  or  the  creation  of  reversible  heat  by  means  of  elec¬ 
tric  current  flow.  The  practical  applications  of  these  phenomena,  such 
as  in  the  cases  of  thermal  refrigeration  and  the  generation  of  electri¬ 
cal  power  from  heat,  has  prompted  a  great  deal  of  research  in  the  last 
decade  in  the  hope  of  optimising  various  thermoelectric  devices.  Theo¬ 
retical  considerations  of  the  results  obtained  from  these  and  more  fun¬ 
damental  investigations  have  revealed  a  rather  complex  interaction  be¬ 
tween  conduction  electrons  and  thermal  currents,  i.e.  electrons  and 
phonons  (1). 

The  purpose  of  this  investigation  is  to  present  a  fundamental 
analysis  of  the  factors  involved  in  the  generation  of  the  Thomson  heat . 
This  purpose  can  be  subdivided  into  two  major  considerations: 

1.  Complete  development  of  the  various  heat  end  electric  flow 
equations  governing  the  evolution  of  heat,  and 

2.  Presentation  of  two  essentially  new  techniques  for  measuring 
the  Thomson  coefficient,  based  on  these  fundamental  equations. 

With  the  application  of  a  tesperature  gradient  to  a  conductor,  in 
the  absence  of  an  electric  field,  not  only  will  the  gradient  generate  a 
flow  of  heat,  but  also  an  actual  sleet ric  current.  In  order  to  observe 
this  current,  however,  it  is  necessary  to  fora  a  complete  circuit  con¬ 
sisting  of  two  dissimilar  metals,  as  the  symmstry  of  a  monometallic 


system  annihilates  any  current  which  is  produced.  The  electric  current 
developed  in  the  circuit  by  the  thermal  gradient  will  depend  upon  the 
resistance  of  the  vires  making  up  the  system.  Alternately  one  could 
avoid  the  need  of  considering  the  resistance  of  the  wire  by  observing 
the  potential  developed  across  the  terminals.  In  this  way  the  observed 
potential  would  be  a  function  only  of  the  thermal  gradient  and  the  par¬ 
ticular  materials  making  up  the  couple.  The  generation  of  this  thermo¬ 
electric  potential  is  known  as  the  Seebeck  effect,  in  honor  of  the  man 
to  whom  the  discovery  is  attributed  (2).  Of  the  three  common  thermo¬ 
electric  effects,  it  is  probably  the  most  familiar  and  widely  used.  A 
term  often  eoplqyed  in  thermoelectric  calculations  is  the  thermoelectric 
power  (S) ,  which  is  simply  the  derivative  of  the  seebeck  potential  (▼) 
with  respect  to  temperature,  i.e.  dV/dT. 

The  second  thermoelectric  effect  to  be  discovered  was  the  Peltier 
heat  («),  which  is  the  heat  developed  in  a  reversible  manner  at  a  Junc¬ 
tion  of  two  dissimilar  metals  when  an  electric  current  is  passed  from 
one  substance  to  the  other  (3).  The  heat  generated  at  such  a  Junction 
should  not  be  confused  with  the  irreversible  Joule  heat,  which  is  al¬ 
ways  evolved  from  a  conductor  carrying  an  electric  current.  The  Peltier 
heat  is  associated  with  the  transfer  of  charge  from  the  interior  of  one 
substance  to  the  interior  of  the  second,  and  is  not,  as  sometimes  as¬ 
serted,  in  any  way  related  to  surface  charges. 

William  Thomson  (h)  proposed  the  existence  of  a  third  thermoelectric 
effect  on  obtaining  erroneous  results  in  applying  the  laws  of  thermo¬ 
dynamics  to  the  first  two  effects.  Thomson  postulated  that  if  an  electric 
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current  is  passing  through  an  individual  conductor  under  an  imposed  tee- 
perature  gradient,  then  there  should  be  a  reversible  heat  developed  per 
second,  over  and  above  the  irreversible  Joule  heat .  This  additional,  or 
Thoason,  heat  is  a  function  of  the  electric  current,  temperature  gradient 
and  the  Thomson  coefficient  (x),  which  is  a  bulk  parameter  of  the  con¬ 
ductor,  the  same  as  resistivity  and  specific  heat. 

Thomson  also  derived  (5)  the  following  relationshipe ,  often  called 
the  "Kelvin  relations",  between  the  three  nnann  thermoelectric 
parameters: 

m  dS 

T 

(1) 

»  ■  TS 

where  8  -  dV/dT 

These  relations  were  first  developed  by  adding  the  irreversible  Joule 
heat  to  the  reversible  Thoason  and  Peltier  heats .  The  addition  of  re¬ 
versible  and  irreversible  heats  is  contrary  to  general  thermodynamic 
considerations  (6),  and,  therefore,  these  equations  were  not  very  satis¬ 
factory  from  a  theoretical  point  of  view.  A  great  deal  of  effort  by 
leading  theoreticians  (6,7)  has  been  made  to  place  these  equations  an  a 
sound  theoretical  basis;  however,  there  has  not  as  yet  been  a  completely 
satisfactory  development  of  these  relations,  with  the  possible  exception 
of  irreversible  thermodynamic  concepts,  where  the  equations  have  been 
developed  with  a  minimum  number  of  assumptions  (8,9,10).  There  can  be 
no  doubt,  however,  of  the  validity  of  these  equations,  .since  they  have 

been  substantiated  time  and  again  by  experimental  investigations. 
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Although  the  Kelvin  relatione  appear  to  offer  the  simplest  experi¬ 
mental  approach  for  determination  of  the  various  thexmoelectrlc  param¬ 
eters,  i.e.  through  measurement  of  the  thermoelectric  power,  a  careful 
analysis  of  the  Thomson  coefficient  data  suggests  that  a  more  direct 
mathod, which  did  not  involve  the  differentiation  of  a  curve,  should  be 
developed.  For  exaaple,  a  number  of  isothermal  Thomson  coefficient 
versus  composition  curves  prepared  from  the  differentiation  of  Seebeck 
data,  presented  by  Salkovits  (11,12)  and  Rudnitskii  (13),  indicated 
that  t  varied  as  a  function  of  composition  in  the  single  phase  region} 
however,  on  crossing  into  the  two  phase  region,  t  often  became  quite 
erratic  with  changing  composition.  If  t  varies  ss  function  of  compo¬ 
sition,  as  demonstrated  in  a  single  phase  solid  solution,  then  in  a 
two  phase  region  in  which  the  conpoeition  of  the  phases  are  constant, 
and  only  the  macro -volumes  are  changing,  one  might  expect  the  varia¬ 
tion  of  t  with  composition  to  behave  somewhat  as  Veg&rd's  law  behaves 
at  a  single  phase  -  two  phase  boundary.  That  is,  variation  to  the 
boundary  in  the  single  phase  region  and  nearly  constant  t  in  the  two 
phase  region.  Such  a  discontinuity  is  not  evidenced  from  the  calcula¬ 
tions  of  t  from  the  Seebeck  data.  In  making  aaasuraments  in  regions 
such  as  this ,  or  in  any  area  where  t  would  be  expected  to  change 
rapidly,  calculations  based  on  the  Kelvin  relations  could  be  subject  to 
error,  since  it  is  vexy  difficult  to  make  precise  slope  measurements  on 
rapidly  changing  curves,  i.e.  dS/dT  (cf.  Equation  1). 

A  simple,  but  fairly  reasonable,  model  of  the  Thomson  heat  may  oe 

developed  from  the  atomic  scale,  if  one  assumes  the  free  electron  model 
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of  *  solid  to  be  an  adequate  first  approximation.  In  such  a  model, 
where  no  interaction  between  the  electrons  and  the  lattice  is  assumed, 
the  thermoelectric  current  can  all  be  attributed  to  the  diffusion  of 
the  electrons  from  the  hot  end  of  the  conductor  to  the  cold  end,  with 
a  consequent  release  or  absorption  of  heat.  Such  a  view  has  led  to 
the  common  terminology  of  calling  the  Thomson  coefficient  the  "specific 
heat"  of  electricity  (l4) ,  since  it  expresses  the  quantity  of  heat  ab¬ 
sorbed  or  evolved  per  unit  of  electric  current  in  passing  through  the 
sample.  Theoretical  calculations  of  the  Thomson  coefficient  based  on 
this  consideration  and  employing  Fermi -Dirac  statistics,  has  led  to 
surprisingly  close  agreement  with  experimental  data  for  a  number  of 
the  simpler  metals  (IS)*  e.g.  the  alkali  metals. 

The  idealised  free  electron  model  is  approached  by  a  few  real 
metals;  however,  there  are  few  cases,  if  any,  where  there  is  abso¬ 
lutely  no  interaction  between  the  free  electrons  and  the  lattice. 

This  being  the  case,  one  would  expect  the  thermal  energy  flowing  down 
the  temperature  gradient,  in  the  form  of  lattice  waves  (phonons)  ,  to 
have  some  effect  on  the  electric  current  in  the  system.  This  has  been 
found  to  be  true  in  some  metals  in  specific  te^e return  ranges.  The 
streaming  phonons  are  believed  to  have  a  tendency  to  sweep,  or  drag, 
the  electrons  from  the  hotter  to  the  colder  end  of  the  conductor.  The 
effect  is  commonly  called  the  "phonon-drag  effect"  and  as  would  be  ex¬ 
pected,  it  is  dependent  upon  the  degree  of  interaction  between  the 
lattice  and  the  electrons .  If  there  is  an  interchange  of  heat  between 

the  electrons  and  the  lattice,  the  "phonon-drag  effect"  must  be  related 
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to  the  ability  of  tha  lattice  to  ahsoxb  ha  at,  i.a.  tha  lattioa  apaeifie 
haat,  aa  illuatratad  in  tha  following  equation  for  tha  contribution  of 
"phonon-drag"  to  tha  Those  on  coefficient  (16) : 


c  ■  lattice  specific  haat  par  atom 
8 

a  ■  probability  of  phonon-electron  collision,  as  eoaopared  to 
other  collisions,  i.a.  phonon -phonon,  phonon-impurity 

8  ■  number  of  conduction  electrons  per  atom 

Aa  tha  temperature  approaches  absolute  sero,  tha  phonon -electron 
collisions  be  cone  dominant  over  other  types  of  phonon  interactions, 
causing  a  to  approach  a  maximum  of  one.  On  tha  other  hand,  the  lattice 
specific  haat  is  known  to  decrease  as  aa  absolute  sero  is  approached. 
At  relatively  high  temperatures  phonon-phonon  interactions  beeoan  domi¬ 
nant,  causing  a  to  approach  sero.  Combining  the  effects  of  the  decay  of 
a  at  high  temperatures  and  the  decmr  of  cg  at  low  temperatures ,  one  would 
expect  the  magnitude  of  the  phonon-drag  contribution  to  reach  a  mast  mum 
at  Intermediate  temperatures,  decmring  on  either  side.  This  was  sub¬ 
stantiated  (17)  experimentally  at  temperatures  between  10  and  50 °K,  where 
thermoelectric  powers  of  aa  high  as  several  mi  111  volts  per  degree  hare 
been  observed.  In  the  metal  systems  investigated  to  date,  the  phonon- 
dreg  contribution  to  the  Thomson  heat  becomes  negligible  in  the  vicinity 
of  room  temperature. 

As  in  the  case  with  electrical  properties,  the  thermoelectric 

effects  are  influenced  by  the  presence  of  impurity  atoam  in  the  lattice. 
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When  the  electrons  tend  to  diffuse  from  one  end  of  the  sample  to  the 
other,  they  will  he  scattered  in  route  by  impurity  atoms.  If  the  scat¬ 
tering  cross  section  is  the  same  for  electrons  with  the  higher  energy 
moving  from  the  hot  end  as  that  for  electrons  with  lower  energy  from  the 
cold  end,  then  the  iaqmrities  would  tend  to  have  no  net  effect  do  the 

thermoelectric  properties.  As  would  be  expected,  however,  the  scat- 

* 

terlng  cross  section  of  an  impurity  atom  is  a  function  of  the  energy  of 
the  electron;  thus  one  finds  that  the  scattering  cross  section.  A, 
would  decrease,  or  the  mesa  free  path,  t,  would  increase,  with  increas¬ 
ing  electron  energy.  If  impurities  are  present  in  a  sample  the  thermo¬ 
electric  equations  must  be  modified  to  take  into  account  the  variation 
of  cross  section  with  electron  energy.  The  effect  of  impurity  atoms 
has  been  studied  at  low  temperatures  (18)  where  the  behavior  was  found 
to  follow  an  equation  of  the  form: 


where 


(1  - 


mMSL  • 

din  E  t 

0 


(3) 


c#Jt  ■  electron  specific  heat 

A(E)  ■  the  energy  dependent  scattering  cross  section 
e  *  charge  on  the  electron 
E  ■  electron  energy 

and  the  derivatives  are  evaluated  at  the  Fermi  level  (c  ) .  The  experi- 

o 

mental  data  agrees  reasonably  well  with  this  equation  at  very  low  tem¬ 
peratures  where  other  processes  are  not  dominant. 
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II.  THEORY 


Since  the  proposal  of  the  Kelvin  relations,  the  theoretical  analysis 
of  the  thermoelectric  properties  of  materials  has  been  presented  from 
tvo  different  standpoints;  irreversible  thermodynamics  as  rigorously  dis¬ 
cussed  by  Onsager  (9,19).  de  Oroot  (10),  Denbigh  (20)  and  Prigoglne  (21), 
or  from  the  kinetic  theoxy  of  electrons  through  Boltsman's  equation,  as 
discussed  by  Wilson  (22)  and  Ziman  (l) ,  to  name  only  a  few.  Since  both 
of  these  approaches  lend  a  degree  of  understanding  to  the  Thomson  coef¬ 
ficient  and  its  relations  to  metals  and  alloys,  a  brief  resume  of  a>ach 
is  presented. 

A.  Irreversible  Therrp^— ' i 

In  the  irreversible  thermodynamic  treatment,  the  creation  and 
transport  of  entropy  plays  a  significant  role.  The  irreversible 
thermodynamic  postulates,  as  reviewed  by  Leech  (23),  can  be  intro¬ 
duced  in  the  form  of  an  equation  linking  the  entropy  flow  (J8)  and 
the  entropy  density  (8^)  by  means  of  the  balance  equation: 

where 

■  electric  and  heat  current  flow  vectors 

■  thermodynamic  "forces"  giving  rise  to  currents 
It  is  further  assumed  that  at  or  near  equilibrium: 

Ji  -  f  Lij 

-8— 


(5) 


la  valid,  and  that  tha  Onsager  ralationa  hold 


(6) 


In  rigorous  treatment  a ,  these  are  deduced  from  fundamental  postu¬ 
lates.  Although  there  is  considerable  choice  in  the  form  of  the 
flow  vectors  or  the  appropriate  conjugate  quantities  X^,  once 
chosen  the  arbitrariness  is  removed. 

The  phenomenological  equations  relating  entropy  flow  density 
J  and  particle  current  density  J  ,  are  then  related  according  to 
Domenicali  (2U)  by: 


-  L 


Vy  - 


VT 


(7) 


-  L  Vy  -  L  VT 


es 


ss 


where  L  ■  L  ,  in  accordance  with  Equation  (6).  However,  a  new 
set  of  phenomenological  equations  are  necessary  to  define  heat 
flow  in  terms  of  J#  and  J#.  following  Callen's  (2$)  work, 
Domenicali  arrived  at: 


J.  -  - -  <£>  («„  -»«jn 

(8) 

V  -  M.q  '  (t>  -  »  V  ” 


vhere  the  M  and  £  coefficients  are  related  by 


v 

v 


6Q,  66 


M«  "  2u  M  +  y  M 
q  eq  e 


chemical  potential 
electrochemical  potential 


If  one  considers  the  steady  electronic  current  J  ,  particles/ 

6 

sec-area,  in  a  cylindrical  rod,  with  no  temperature  gradient,  then 
the  ratio  of  Equations  7A  :  7B  reduces  to: 


VJ  'VT»0 
6 


~  S  8* 


(9) 


which  defines  the  "transport  entropy  per  particle".  The  ratio 
of  Equations  8A  :  8B  yields: 


(i) 

VJ  'VT-0 
6 


M  . 

-21  r  Q* 

Mee 


(10) 


which  defines  the  "transport  heat  per  particle".  It  can  be 

2 

further  shown  that  L  ■  o/e  ,  where  c  is  the  isothermal 

66 

conductivity . 

These  terms,  developed  from  fundamental  concepts,  can  be 
directly  applied  in  the  determination  of  the  values  of  the  various 
thermoelectric  coefficients.  Values  of  the  Thomson  coefficient, 
for  example ,  can  be  determined  by  applying  these  concepts  to  the 
equation  developed  to  calculate  the  additional  rate  of  heat  flux 
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required  to  maintain  a  given  temperature  distribution  once  the  elec¬ 
trical  current  ia  turned  off.  The  solution  takes  the  form: 

eAE  •  -  /y  (T  ||r)  J.VTdv  -  e  fy  (^)  At  (11) 

where 

1  38 

t  ■  -  (— )  T  ■  Thomson  Coefficient 

J  ■  electrical  current  density  =  -  e  J# 

AE  ■  the  rate  of  heat  transferred  into  volume  element 
v  to  maintain  temperature  gradient 

In  order  to  determine  x ,  it  is  necessary  to  remove  it  from 
under  the  integral  alga,  which  can  moot  easily  be  accomplished  by 
considering  a  long  homoganeous  section  of  thin  wire.  In  the  steady 
state  situation,  Equation  (11),  then  reduces  to  the  experimentally 
usable  relation: 

AE  •  A  *  It  AT  (12) 

The  reversible  Thomson  heat  (xl  AT)  adds  or  subtracts  to  the  irre¬ 
versible  Joule  heat  (I^R),  depending  upon  the  direction  of  electric 
current  flow  with  respect  to  the  te^erature  gradient.  The  sign 
convention  followed  in  this  investigation  is  that  established  by 
MacDonald  (17)*  in  which  the  positive  sign  indicates  the  electrical 
current  is  flowing  down  the  thermal  gradient,  i.e.  from  Tg  to 
where  Tg  >  T^,  and  the  negative  sign  indicates  the  electrical 
current  is  flowing  up  the  thermal  gradient. 

Equation  (12)  gives  a  very  satisfactory  derivation  from  funda¬ 
mental  considerations  in  the  steady  state,  which  can  be  used  in 


direct  measurement* 


B.  Kinetic  Theory 

The  kinetic  theory  approach  to  trensport  properties  in  solids 
is  introduced  through  the  assumption  that  the  electron  distribution 
in  the  free  electron  gas  nay  be  represented  by  sn  equation  de¬ 
veloped  by  7.1  nan  (1) ,  which  is  of  the  form: 

du  »  f(Vvy,vt  ;  x,y,a  ;  t)dvxdvydv8dxdyds  (13) 

where  du  is  the  nusfcer  of  electrons  at  tine  t_  in  the  volume  elenent 
dx,  dy,  da,  with  velocities  between  and  (vx  +  dvx).  Accepting 
the  most  elementary  cose  and  allowing  Equation  (13)  to  define  "f", 
the  distribution  function  of  the  electrons,  then: 

/{/  f  dvxdvy drx  •  Kxjr.x  ;  t)  (lU) 

where  H  is  the  density  of  free  electrons  at  tine  The  distri¬ 
bution  function,  f,  has  the  equilibrium  value,  fQ,  in  pn  isothermal 
conductor  rod,  not  under  the  influence  of  electric  or  nagaetlc 
fields,  and  is  given  by  the  Fermi -Dirac  function: 

f o  ■  2  (jj)3  (eE  ”  *  ♦  l)’1  (15) 

in  which  E  is  the  energy  of  on  electron  with  velocity  v  ,  v  ,  v 

A  v  ” 

(often  presented  as  k  vectors),  snd  e  is  the  chemLcsl  potential 
of  the  electrons. 

The  Boltsmann  equation  is  presented  formally  as: 

_  /££)  B  (££)  (i<\ 

'dt' fields  'dt ^collisions  lAO; 
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where  the  left  hand  aide  represents  the  effect  of  external  flelda 
on  the  electron  without  interference*  and  the  right  hand  side 
represents  the  effects  of  electron  collisions  on  the  distribution 
function. 

For  the  simplest  case  where  the  electron  flow  is  due  to  an 
electric  field  and  a  temperature  gradient  dX/dx*  and  all  collisions 
tend  to  restore  equilibrium  at  a  linear  rate  of  restoration*  which 
is  proportional  to  the  departure  from  fQ*  Equation  (16)  nay  then 
be  written: 


3v 


!Ss.  ffiT  .  *' 

"  41  4T  * 


f-f 


(17) 


ratures*  T  >  8^,  or  when  scattering  is  due  to 


where  it  is  assunad  that  t  (c)  is  the  constant  restoration  rate* 
irrespective  of  the  type  of  deviation  from  equilibrium.  The 
assumptions  leading  to  the  restoration  rate  are  approximately  true 
only  at  high  ti 
inpurity  atone. 

Utilising  these  assumptions*  Leech  (23)  derives  a  set  of  phe¬ 
nomenological  equations  based  on  kinetic  theory  which  are  siailar 
to  the  irreversible  themodynanic  phenomenological  Equations  (7) 
and  (8): 


j  .  gglT.Ill  fi  -•,'(«)  d3r£. 

x  n  ^x  T  8T  '  Tx  *o  dx 


(18) 


Q  .  .  £1— LlL  /  T  v2  ffo  .3  a  _  SLlii  /  v  f  t  d3v  ® 

Hx  2  1  V  9r  a  T**x  2  J  V  ro  a  T  dx 


(19) 
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vhsre 


H  ■  /  f  d3v 


electron  density 


Relating  these  equations  to  those  developed  from  irreversible  therao- 
dynsad.es.  Leech  redefined  the  Those  on  coefficient  in  terse  of  the 
kinetic  theory  es: 


T  L  (i  Isl  .  i£} 
1  8T  ^3  e  8T; 


(20) 


vhere: 

C0l  ■  electronic  heat  capacity 
g  ■  cheaical  potential/unit  conduction  charge 
For  metallic  conductors,  g  is  the  Fermi  energy  per  electron  per 
electron  charge,  and  if  the  following  free  electron  model  re¬ 
lations  hold: 


*  •  vk^T2 

«  ■  ~ nr 


12  c. 


(21) 


vhere 


c 


e 

o 


£  <^2/3 


Feral  level  act  0°K 


and 


C™  *  w^k^T^e*  ■  electronic  heat  capacity 

then  the  Thomson  coefficient  is  related  directly  to  the  electron 
heat  capacity  per  electron: 


T 


(22) 


which  vas  originally  suggested  by  Thomson 


C.  Heat  Balance  Equations 

The  aquations  developed  in  the  previous  sections  establish  e 
means  of  date  mining  the  Thomson  coefficient  from  first  principles. 
The  resulting  equations,  however,  are  derived  from  equilibrium  con¬ 
siderations  in  a  specific  conductor}  and  as  such,  do  not  give  a 
complete  analysis  of  all  the  heats  involved  in  the  system. 

The  general  equation  vhich  is  employed  for  problems  vhich  in¬ 
volve  heat  flows  into,  and  out  of  a  system  is: 

-  Qa  ♦  ^  (23) 


where: 

Qj  ■  rate  of  heat  input  to  the  system 
Qa  *  rate  of  heat  absorbed  by  the  sample 
■  rate  of  heat  leaving  a  unit  volume. 

One  should  note  that  these  are  rate  of  heat  transport,  and  not, 
as  often  stated,  bulk  quantities  of  heat.  To  gain  a  complete 
understanding  of  the  phenomena  taking  place  within  a  metal  rod 
subjected  to  the  various  energy  flows,  it  is  necessary  to  under¬ 
stand  how  a  given  rate  of  heat  input  can  give  rise  to  one  equi¬ 
librium  tasperature  and  a  different  rate  of  input,  give  rise  to 
a  different  equilibrium  te^erature.  ibis  is  briefly  described 
in  Appendix  A. 


Seits  (26)  has  given  a  detailed  study  of  the  Th< 
efficient  using  Equation  (23)  in  the  form: 

£aIi 

dt  Ix  x  3x 


co- 


(24) 
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Wher*  dH/dt  is  the  rats  at  which  hsat  accumulates  in  a  unit  volume 

of  conductor  carrying  both  electrical  and  thermal  currants.  Ths 

ten  (I  E  )  in  this  equation  is  related  to  the  electrical  work  done* 
3Cx 

and  the  term  (— “)  is  related  to  the  rate  of  heat  flowing  out  of  a 
unit  volume .  The  equation  can  be  written  in  a  more  familiar  font 
by  writing  E^  and  Cx  as 


E 


..L/iL!  _  i  s_  (tl)  a.  _  !i  i_  a. 

#2r2K1Ax  e  dT  '  T'  dx  ^  T»  dx 

c  ,.i  ij  ,  i  /r [*3ici-|‘z)  « 

Sc  ^eAx  3’»'  TK^  1 


dx 


where: 


3/2  #  \/2  ;o  fl  1  to  ** 


(25) 


(26) 


nQ  *  nuaber  of  electrons  per  unit  volume 
• 

c  ■  energy  of  Fermi  electrons  at  0°K 
o 

t  ■  Man  free  path  of  electrons 
f  ■  distribution  function 

substituting  Equations  (25) ,  (26)  into  Equation  (24),  one  obtains 


«/*t- !*,♦?*  *  -  V  *  I?  <«  S> 


(27) 


The  first  term  to  the  right  of  the  equality  in  Equation  (27) 
represents  the  Joule  heat,  and  the  third  is  due  to  that  arising 
from  the  conduction  of  heat  in  the  material.  Both  of  these  tei 
are  independent  of  the  relative  direction  of  the  electric  and 


thermal  current* ,  while  the  second  tens  represents  a  rate  of  accu¬ 
mulation  of  heat*  dependent  upon  the  relative  direction  of  the  two 
flows. 


Rewriting  the  second  term  of  Equation  (27),  heat  will  be  pro¬ 


duced  at  a  rate: 


rS)  .  fs.  4_(Jk.ill  i  a. 

^ dt  t  e  dl  VTKX  T'  A  dx 


which  can  now  be  recognised  as  the  Thomson  heat,  where  the  negative 
of  the  bracketed  term,  namely 


.1  4_  ( Ji  .  si)  .  t 

e  dET  'DL  T'  T 

is  equivalent  to  the  Thomson  coefficient.  This  equation  can  be 
further  refined  leading  to  Equation  (29): 


•  3  k  »(.  )J 


where  the  quantity  1  (eo)/l(eQ)  can  be  roughly  evaluated  using  a 
theory  of  the  mean  free  path. 

Simplifying  the  nomenclature ,  Equation  (27)  can  be  rewritten: 

i.£*k$  <»: 

Under  steady  state  conditions,  i.e.  dB/dt  ■  0,  this  becomes 
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This  differential  equation  it  solved  in  Appendix  B.  The  resulting 
expressions  ere  then  used  in  s  direct  determination  of  the  Those  on 
coefficient,  for  both  general  and  special  cases. 


III.  OTHER  EXPERIMENTAL  WORK 


The  method  most  commonly  used  for  the  direct  measurement  of  the 
Thomson  coefficient  is  based  on  a  method  perfected  by  LeRoux  (27) ,  and 
used  with  slight  modifications  by  Berg  (28),  Hettleton  (29),  Young  (30), 
Smith  (31),  Borelius  (32),  and  Systrom  (33).  The  actual  experimental 
arrangements  used  to  establish  the  tesperature  gradient  and  design  of 
the  apparatus  varies  slightly  from  one  investigation  to  another,  the  main 
differences  lying  in  the  means  of  detecting  the  temperature  changes  and 
gradients  in  the  sample. 

According  to  Hystrom  (33).  on  subjecting  a  homogeneous  sample  to 
a  constant  thermal  gradient  and  electrical  current,  the  governing 
equation  can  be  written: 

heat  produced  »  J2 p  ‘  Jt  ^  (32) 

If  the  direction  of  current  is  reversed,  the  first  term  denoting  Joule 
heat  will  not  change,  but  the  second,  denoting  the  Thomson  heat,  will 
change  sign.  If  the  current  strength  is  altered  by  a  suitable  amount, 

AJ,  in  connection  with  its  commutation,  the  total  effect  will  be  un¬ 
altered.  By  equating  the  two  equations,  one  can  calculate  x. 

The  mathematical  development  of  the  foregoing  begins  with  the  cur¬ 
rent  flow  in  a  direction  with  respect  to  the  thermal  gradient  such  that 
the  Thomson  heat  is  added  to  the  Joule  heat,  l.e. 

A  ♦  *  £ 

then,  with  a  reversal  and  subsequent  change  of  current,  one  can  return 
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to  the  original  equilibrium,  i.e.i 


A  ♦  jt  g-  (J  +  AJ)2p  -  (J  ♦  AJ)t  — ■  (33) 

expanding  and  rearranging: 

2tJ  AJx  2JAJp  ♦  dA  (3U) 

and 


Thus  t  can  be  calculated  directly  from  the  experimental  conditions. 

Early  experiments  employed  a  bolometer  coll  to  detect  the  tempera¬ 
tures  of  the  sample.  In  this  method  a  copper  resistance  thermometer 
coll  vas  wound  around  the  test  specimen  as  shown  in  Figure  1,  and  the 
change  in  resistance  is  detected  by  a  very  sensitive  bridge  circuit . 
Since  the  variation  of  resistance  of  copper  with  temperature  is  known 
quite  accurately,  the  temperature  of  the  simple  at  the  region  in  ques¬ 
tion  is  readily  determined.  At  thermal  equilibrium,  as  detected  by  the 
bolometer  coil,  the  current  vas  reversed  and  then  varied  until  the  tem¬ 
perature  of  the  sample  vas  Identical  to  that  established  vith  the  orig¬ 
inal  current  flow. 

More  recent  investigations  (13) »  have  used  a  similar  analysis, 
only  esplcying  thermocouples  as  the  sensing  elements.  The  thermocouples 
were  either  taped  directly  to  the  side  of  the  specimen  or  inserted  into 
small  surface  holes  drilled  into  the  specimen. 

Hy stroa  (33)  has  used  a  microthermopile  to  measure  both  the  tem¬ 
perature  gradient  and  the  temperature  changes  associated  vith  reversal 
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of  currant  flow.  The  projector  and  the  thermopile  were  built  such  that 
the  entire  "microscope"  could  be  moved  to  obtain  a  maximum  emf  out  of  the 
projected  image  of  the  saaple.  The  output  of  the  thermopile  vae  meaaured 
with  a  very  accurate  potentiometer.  Again  the  procedure  for  obtaining 
the  Thomson  coefficient  was  the  aame  as  that  described  previously. 

This  basic  technique  is  quite  valid  and  has  provided  a  great  deal 
of  basic  data  on  the  Thomson  coefficient.  Analyses  have  been  made,  how¬ 
ever,  in  which  Equation  (32)  has  been  referred  to  as  a  certain  amount  of 
heat,  which  causes  a  change  in  temperature  of  the  conductor.  Techni¬ 
cally,  this  is  incorrect,  since  it  is  not  a  quantity  of  heat  at  all,  but 
the  rate  of  heat  input  to  the  entire  system,  none  of  which  is  absorbed 
by  the  conductor  at  equilibrium,  which  although  not  frequently  mentioned, 
is  the  actual  case  considered  in  this  type  of  analysis.  If  one  does 
carry  out  this  technique  it  is  also  found  that  dT/dx  is  not  simply  the 
difference  in  temperature  between  the  end  sinks  divided  by  the  length  of 
the  sample,  as  is  often  inferred,  but  varies  throughout  the  sample  and 
also  changes  value  with  reversal  of  current.  The  only  position  at  which 
dl/dx,  for  current  flow  in  one  direction,  is  the  same  as  that  for  current 
flow  in  the  opposite,  is  exactly  at  the  midpoint  of  the  sample.  A  com¬ 
plete  analysis  of  this  system  is  given  in  Appendix  B. 

In  general  then,  this  technique  is  quite  valid  for  obtaining  values 
of  the  Thomson  coefficient.  Some  of  the  current  literature  (13,33),  how¬ 
ever,  gives  an  incomplete  development  of  the  theory  governing  the  analysis, 
which  can  make  it  quite  difficult  to  obtain  a  basic  understanding  of  the 
factors  involved  in  calculating  the  Thomson  coefficient. 
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Lander  ( 31*)  has  used  a  rather  unique  approach  for  a  direct  Measure¬ 
ment  of  the  Thomson  coefficient.  He  assumed  that  the  conduction  heat  out 
of  the  end  of  a  wire.  Ignoring  radiation  losses,  is  equal  to  the  Joule 
heat  plus  the  Thomson  heat  developed  between  that  end  and  the  surface  of 
maximum  temperature.  A  sketch  of  his  apparatus  is  shown  in  Figure  2. 

The  two  ends  of  the  sample  are  maintained  at  the  same  temperature  and  the 
center  of  the  sample  is  heated  electrically.  It  will  be  noted  that  with 
an  electrical  current  flow  (I)  in  one  given  direction,  the  current  flow 
in  one-half  of  the  sample  will  be  in  the  same  direction  as  the  thermal 
gradient  and  in  the  other  half,  against  the  gradient. 

In  the  absence  of  Thoomon  heat  the  conduction  heats  flowing  out 
the  ends  of  an  Isotropic  wire  are  equal  if  these  ends  are  at  the  earns 
temperature.  Consider  a  wire  with  both  ends  held  at  absolute  temperature 
Tq  and  heated  electrically  to  a  marl  mum  tasperature  T*.  Let  the  elec¬ 
trical  potential  at  any  point  be  indicated  by  +  and  the  potential  at  the 


two  ends  be  0  and  8.  The  heat  flowing  per  second  through  the  surface 
with  0  potential  and  temperature  Tq  is,  in  the  steady  state: 


*1  "  ♦  I  ♦  I  i\ 


TdT 


(36) 


*  .  • 

where  +  is  the  potential  at  T  and  I  is  the  electrical  current  flowing 
through  the  wire.  Through  the  surface  at  potential  8  and  temperature 


Tq  the  heat  flowing  per  second  is: 


S 


(8  -  ♦*)  I  -  I  /J  TdT 
o 


(37) 
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To  Power  Circuit 


FIGURE  1  -  Apparatus  first  usad  to  measure  Th onion  Coafflciant 

Current  Leads  .  . 


FIGURE  2  -  Apparatus  usad  by  Lander  to  measure  Thoaaoo  Coafflciant 
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The  difference  between  these  divided  by  the  total  heat  output  is 


♦<|)/*\«  <*» 

Through  mathematical  considerations  end  certain  simplifying  assump¬ 
tions  this  equation  finally  reduces  to: 

(I»a  4*  *<***-  ^)1/s  o  (») 

o 

which  nay  be  s imply  evaluated  if  the  Thoason  coefficient  is  proportional 
to  the  temperature.  Unfortunately,  this  assumption  is  not  generally 
Justified,  as  indicated  by  MacDonald  (35).  Tor  evaluation  then,  it  is 
necessary  to  represent  either  the  right  or  left  side  of  the  equation  by 
a  power  series.  The  usual  amthod  is  to  represent  the  Thoason  coefficient 
by  a  power  series: 

x  ■  a  ♦  0T  ♦  yT*  ♦  0T3  ♦  .... 

It  is  then  found  that 

(^q)  (§)  ■  ««**  ♦  wr*2  ♦  yct*3  ♦  ear*1*  ♦  ....  (to) 

e 

where  T  is  aatheantlcally  evaluated  assuming  the  validity  of  the 
ttLedaaann-Frans  Law.  The  terns,  a,  0,  y,  and  0  are  evaluated  fron  the 
experimental  data  of  ^  ,  and  the  Thoason  coefficient  is  calculated  in 
turn  from  these  values. 

A  great  deal  of  experimental  work  has  been  performed  by  many 

different  investigators  (36,  37)  on  indirect  means  of  determining  the 

Thoason  coefficient.  All  of  the  techniques  e^loy  the  Kelvin  Relations 

to  calculate  the  various  thermoelectric  parameters  by  measuring  the 
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change  of  Saabeck  potential  ( V)  with  respect  to  te^erature  (thermoelec¬ 
tric  power)  and  then  determining  the  Thoaaon  coefficient  with  Iquation 
(1). 

It  la  impossible  to  describe  all  the  significant  contributiona 
made  to  the  knowledge  of  the  race  le  Ctrl  ca  through  these  maasureaants ; 
however,  one  group  which  has  contributed  significantly,  especially  at 
low  temperatures ,  was  that  headed  by  the  late  D.K.C.  MacDonald  (IT)* 

In  making  measurements  at  very  low  temperature,  MacDonald's  group 
has  mads  great  use  of  the  fact  that  certain  elements  and  compounds  are 
super-conducting  to  relatively  high  teaperatures.  Since  in  the  super¬ 
conducting  temperature  range,  these  materials  exhibit  aero  thermo¬ 
electric  properties ,  any  Beebe ck  voltage  measured  between  this  and 
another  material  will  all  be  due  to  the  second  material.  In  this  way 
absolute  Seebeck  potential  can  be  measured  directly. 

The  low  temperature  apparatus  used  by  Christies,  Jan,  Pearson, 
and  Templeton  (38)  for  the  determination  of  the  absolute  scale  of  lead, 
is  similar  to  that  used  in  low  temperature  thermal  conductivity  mea¬ 
surements  (39)  •  In  their  approach  a  couple  was  made  between  the  lead 
and  the  superconducting  Ib^8n,  end  the  temperature  gradient  Imposed  by 
electrical  means.  The  Seebeck  voltage  was  measured  quite  accurately 
with  a  very  sensitive  potentiometer  circuit.  The  temperature  was  then 
varied  and  a  plot  of  thermoelectric  power  vereua  te^erature  was  marie, 
from  the  slope  of  this  curve,  the  Thomson, coefficient  could  be  calcu¬ 
lated  from  ths  Kelvin  Relations . 

The  results  obtained  at  low  teaperatures  on  various  materials  are 

-25- 


very  interesting  and  rather  surprising  (1*0) ,  in  that  the  thermoelectric 
power  varies  a  great  deal,  even  to  the  extent  of  changing  sifi.  Hie 
results  do  not  agree  with  the  classical  theory  of  themoelectrios  (IT)  • 
and  as  such  have  provoked  a  great  deal  of  diacussion  (1,17, bl)  as  to  the 
nature  of  interactions  taking  place  at  low  temperatures . 


IV.  ZXPSRIMBrTAL 


Ab  previously  indicated,  the  experimental  apparatus  necessary  for 
the  direct  msasuremant  of  the  Thoaeon  heat  oust  include  mbi  for 
applying  a  thermal  and  a  D-C  potential  gradient  ,  and  adequate  thermome¬ 
try  to  meaaure  both  the  thermal  gradient  and  the  thermal  changes  asso- 
ciated  with  the  reversal  of  currant  flow.  In  order  to  avoid  a  long 
tedious  analysis  of  heat  losaes(  in  which  certain  simplifying  assump¬ 
tions  sust  be  aade,  the  system  should  also  incorporate  means  of  reducing 
radial  heat  losses  from  the  sample  to  the  surroundings.  A  calorimeter* 
to  operate  in  the  teaperatura  range  of  *  10C°C,  vas  designed  and  con¬ 
structed  to  optimise  each  of  these  requirements. 

Also  incorporated  in  the  experimental  considerations  vas  a  ■  angle 
design  intended  to  improve  on  previous  methods  of  detecting  the  thermsl 
changes  associated  vlth  the  Thomson  heat.  The  standard  techniques 
generally  employed*  1  .a .  thermocouples  velded  to  the  sides  of  the  speci¬ 
men,  or  the  standard  optical  techniques,  do  not  give  valid  thermal  sea- 
aureaents.  The  use  of  thermocouples  in  this  manner,  for  example, 
chengse  the  temperature  distribution  in  the  area  of  oontact,  while  op¬ 
tical  measurements  at  low  temperatures  are  subject  to  errors  introduced 
by  diffraction  and  absorption.  Tor  this  reason  the  sables  in  this  in¬ 
vestigation  were  designed  with  a  vary  small  axial  hole  extending  the 
full  length  of  the  specimen,  such  that  a  veil  insulated  thezmooomple 
could  be  essentially  embedded  in  the  samples. 
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A.  Calorimeter 

The  outer  shell  (F)  of  the  calorimeter  wee  constructed  fro*  e 
3"  cylinder  of  solid  copper  as  shown  in  Figure  3.  la  order  to  pre¬ 
vent  the  current  froa  flowing  through  the  ealorlaeter  shell,  rather 
than  through  the  specimen,  the  two  end  sections  of  the  ealorlaeter 
were  electrically  Insulated  froa  each  other  by  constructing  the 
ealorlaeter  in  two  sections,  separated  by  a  heavy  rubber  gasket  (B), 
and  held  together  by  nylon  bolts  (M) . 

Massive  copper  cylinders ,  designed  to  slip  fit  into  the  end 
shells  of  the  ealorlaeter,  served  as  temperature  control  sinks  (C) . 
Kadi  sink  contained  twelve  evenly  spaced,  1/U"  holes,  approximately 
1/8"  froa  the  outside  circumference,  extending  the  length  of  the 
cylinder  (l),  into  each  of  which  were  inserted  hollow  alumina  in¬ 
sulating  tubes.  Besting  coils  of  0.030"  chrcssel  heating  wire  were 
in  tun  placed  in  the  Insulators .  The  twelve  helical  coils  were 
connected  in  series  and  existed  the  ealorlaeter  through  Insulated 
seals  (D)  in  the  top  of  the  ealorlaeter  to  a  0  to  100  volt  varlac. 

A  temperature  control  thermocouple  well  (I)  was  placed  near  the 
inside  edge  of  each  cylinder  to  insure  marl  nun  temperature  control 
of  the  spedaen.  The  chroael-eluael  thermocouple  in  each  of  the 
sinks  existed  through  insulated  seals  to  the  temperature  control" 
ling  circuit.  The  copper  top  plates,  to  which  the  insulated  vacuum 
seals  were  soldered,  were  bolted  to  the  hodr  of  the  ealorlaeter. 

The  teaperature  control  sinks  were  designed  to  hold  a  1" 

QD  x  1"  cylindrical  saaple  heat  sink,  as  shown  in  Figure  k.  The 
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FIGURE  3  -  Calorimeter  for  Thonson  Coefficient  determination 
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sinks  themselves  ware  mads  of  the  material  under  teat  in  order  to 
minimise  any  Peltier  effects  between  the  sample  rod  and  the  heat 
sink.  A  hole,  the  else  of  the  test  sa^le ,  was  bored  in  each  of 
the  sample  heat  sinks  to  allow  maximum  thermal  contact  with  the 
specimen. 

The  test  saaples  were  1U9.I  os  long,  with  an  outside  diameter 
of  approximately  1*.95  mm.  A  thermocouple  hole,  varying  from  0.5 
to  1.0  mm  in  diameter,  depending  upon  the  type  of  specimen  being 
tested,  was  either  milled  or  formed  in  the  specimen  during  casting. 
An  electrically  insulated  thermocouple  was  inserted  in  the  hois  to 
measure  the  temperature  gradient  and  the  change  in  temperature 
associated  with  the  change  in  currant. 

During  operation,  the  surface  of  the  specimen  was  maintained 
under  a  vacuum  to  minimise  radial  heat  losses  from  the  sample  to 
the  surroundings.  This  necessitated  a  number  of  "O"  ring  seals  as 
illustrated  in  Figure  3  and  k.  It  was  also  necessary  to  use  a 
vacuum  sealing  varnish  (Qlyptol-General  Electric)  around  the  ends 
of  the  saqple  to  prevent  vacuum  leaks  along  the  interface  between 
the  sample  and  the  saaple  heat  sinks.  To  further  minimise  the  heat 
losses  from  the  sample  by  radiation,  a  silvered  glass  tube  (I)  was 
placed  around  the  specimen. 

B.  Thermocouples 

The  preliminary  Investigations  using  chromel-aluaml  thermo¬ 
couples  for  temperature  measurements  of  the  sample  yielded  output 

emfs  which  were  quite  erratic.  Two  sources  of  instability  were 
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determined,  i.e.  inhomogeneitiee  in  the  thermocouple  vires  and  stray 
electromagnetic  fields  around  the  exposed  thermocouple  leads. 

Inhomogeneitiee  in  the  chronel  and  aluael  thermocouple  vires 
vers  found  to  give  rise  to  stray  Seebeck  voltages  vhen  the  vires 
were  exposed  to  snail  temperature  variations  along  its  length.  The 
inhonogeneitiea  theme  elves  act  as  point  sources  of  voltages,  vhich 
could  be  detected  by  the  very  sensitive  measuring  circuit .  If,  for 
example*  the  lead  vires  were  subjected  to  stray  air  currents,  this 
vas  sufficient  to  give  rise  to  temperature  variations  along  the  vire, 
vhich  could  be  detected  by  the  sensing  circuit .  The  addition  of 
thermal  insulation  to  the  vire  tended  to  reduce  this  effect,  but  not 
of  sufficient  magnitude  to  give  satisfactory  measurements. 

Hultgren  (42)  Indicated  that  selected,  annealed  copper-cons tan- 
tan  thermocouples  vere  of  sufficient  homogeneity  to  allow  stable  mea¬ 
surements  of  very  small  temperature  differences.  The  actual  vires 
to  be  used  for  the  thermocouple  vere  carefully  selected  by  exposing 
small  sections  of  the  vire  to  sharp  thermal  gradients  while  the  two 
ends  of  the  vire  vere  connected  to  a  Nanovoltmeter  (scale  to  10”^ 
volts).  Those  lengths  which  gave  the  minimum  output  voltage  under 
the  thermal  gradient  vere  finally  selected  for  the  actual  thermo¬ 
couples  . 

In  order  to  obtain  the  smallest  diameter  hole  possible  in  the 

saaple,  the  thermocouple  bead  vas  formed  as  a  butt  veld,  such  that 

the  two  leads  could  exit  the  sample  in  opposite  directions,  as  shown 

In  Figure  5.  The  bead  itself  vas  cars  fully  filed  to  approximately 
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FIGURE  5  -  Tcapcratur*  recording  syatca 
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the  ease  size  as  the  vire  end  electrically  Insulated  fro*  the  sample 
by  forcing  thin -vailed  teflon  shielding  over  the  thermocouple  vire. 
Since  high  D-C  currents  flowed  axially  around  the  thermocouple,  ex¬ 
treme  care  vas  exercised  in  assuring  that  no  extraneous  currents  were 
transmitted  to  the  thermocouple  circuit.  This  vas  tested  by  placing 
the  thermocouple  in  the  sample  and  impressing  high  currents  through 
the  sample  and  examining  the  output  for  variations  that  could  be 
attributed  to  the  changes  in  the  current .  If  any  changes  of  this 
nature  were  detected,  the  thermocouple  vas  removed  and  the  insula¬ 
tion  replaced. 

Since  variations  in  the  thermocouple  emf  of  the  order  of  0.0$ 
microvolts  were  to  be  investigated,  precautions  which  are  not  nor¬ 
mally  necessary  in  standard  thermal  measurements  were  required.  TO 
prevent  interference  from  the  ever  preeent  electromagnetic  fields 
generated  by  control  devices ,  all  instrument  leads ,  as  well  ss  the 
exposed  section  of  the  thermocouple,  required  very  careful  shielding 
to  ground.  In  detecting  the  small  differences  in  temperature  asso¬ 
ciated  with  the  Thomson  heat,  the  cold  Junction  also  becomes  more 
significant.  Small  voltage  fluctuations  associated  with  the  melting 
ice  were  at  first  detected  in  the  measurements.  Substitution  of  a 
deep  devar  flask,  packed  very  tightly  with  crushed,  distilled  ice 
and  a  minimum  of  distilled  water,  provided  a  very  stable  cold  Junc¬ 
tion.  To  further  insure  against  any  possible  short  circuits  by 
electrochemical  means  between  the  two  thermocouples ,  both  cold  Junc¬ 
tion  terminals  were  electrically  insulated. 
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C.  Electrical  Circuit  and  Controls 

The  output  eaf  from  the  thermocouple  vas  counterbalanced  by 
the  output  eaf  of  a  Lin  deck  microvolt  source  (Model  IS) ,  such  that 
all  or  any  part  of  the  thermocouple  signal  could  be  cancelled,  as 
shown  in  Figure  5.  In  this  Banner,  any  eaf  variation  over  a  de¬ 
sired  temperature  range  acted  as  a  source  for  a  Kelthley  Nano¬ 
volt  meter  (Model  148),  such  that  differential  changes  in  the  order 
of  3  x  10“5oC  could  be  detected. 

The  one  volt,  full  scale  output  from  the  H ano volt ae ter,  re- 

D 

gardless  of  the  scale  (1,  3,  10,  300  ...  x  10  volts),  vas  reduced 
by  a  fixed  voltage  divider  network,  the  output  of  which  acted  as 
a  source  for  a  Brown  Electronic  Recorder  (Model  Ho.  13301876).  The 
recorder  was  modified  according  to  Kanda  at  al  (43)  to  allow  for  a 
variable  range  froai  0  to  50  av  and  a  variable  aero  point  adjustment 
frost  o  to  20  av.  In  operation  the  zero  point  and  the  span  of  tba 
Brown  Recorder  were  calibrated  vith  one  scale  of  the  Hanovoltaster . 
Since  the  full  scale  output  of  the  Hanovoltaster  was  constant,  re¬ 
gardless  of  the  scale,  macro  or  micro  variations  In  temperature 
could  be  recorded  directly  on  the  recorder. 

D.  Power  Supply  to  S«mie 

The  electrical  system  used  to  supply  the  current  for  the 
specimen  is  shown  in  detail  in  Figure  5.  The  currant  was  supplied 
by  three  fully  charged  4  volt  cadaiua  batteries  (Cadaidyne)  in 
parallel  vith  an  electro  source  D-C  power  supply  (Nodal  EFB).  The 

current ,  regulated  by  a  4  ohm,  1000  watt  Ohaite  power  rheostat 
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(Motel  U) ,  flowed  through  the  calibrated  reeietor  of  90  sill  aangan 
wire,  through  the  current  reverting  twitch  and  the  sample.  Contact 
resistance  at  each  end  of  the  terminals  was  naintained  act  a  mini  sue. 
The  12  gauge  copper  wire  used  in  the  circuit  was  connected  such 
that  the  circuit  resistance  was  independent  of  the  direction  of 
current  flow.  The  current  in  the  circuit  was  constantly  Monitored 
during  a  particular  investigation  by  using  the  voltage  drop  across 
the  known  aangan  an  resistance  as  an  enf  source  for  a  second  Modified 
Brown  Electronic  Recorder.  The  recorder  ves  standardised  with  the 
output  from  the  Linda ck  microvolt  source.  Current  changes  in  the 
saaple  of  less  than  .01  «ps  could  be  detected  with  this  circuit. 

The  ewe  ter  (Meetinghouse  -  Model  PX-h)  wee  used  to  roughly  select 
the  current  to  be  tested. 

Various  current  reversing  switches  were  tested  in  the  circuit j 
however,  it  was  found  that  a  a  lapis  heavy  duty  double -pole -double* 
throw  knife  switch.  Modified  such  that  the  poles  were  physically 
within  a  few  millimeters  of  each  other  at  the  open  circuit  position, 
gave  practically  aero  switching  tine  and  a  minimum  of  switching 
noise.  In  order  to  insure  that  the  current  did  not  change  value  in 
the  switching  process,  identical  short  pieces  of  wire  were  care¬ 
fully  silver  soldered  into  the  circuit,  so  that  the  resistance  of 
both  paths  were  Identical.  It  was  found  that  the  currant  changes 
were  almost  sero  during  the  reversal  of  currant. 

x.  SfegflMdL -Sadias* 

The  temperature  gradient  in  the  sample  was  established  by  the 
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two  end-tenpe rature  control  hoot  oinks  of  the  calorimeter  as  des¬ 
cribed  in  the  calorimeter  section.  The  temperature  of  the  sinks 
vers  originally  controlled  by  a  multipoint  Tsbilaug  Celect-Ray 
Controller  (Model  U8333) *  as  activated  by  the  respective  control 
thermocouples  placed  in  the  end  heater  sinks.  Various  temperature 
gradients  for  the  sample  could  be  chosen  by  presetting  the  de¬ 
sired  temperature  on  the  Celect-Ray  unit.  In  the  first  series  of 
experiments  it  vas  noted  that  the  sample  thermocouple  vas  detect¬ 
ing  slight  changes  in  temperature ,  corresponding  to  the  on-off 
cycle  of  the  heaters .  For  this  reason  the  Thomsen  heat  data  to 
follow  was  obtained  with  the  heater  sinks  uader  static  heating 
conditions;  that  is,  the  high  temperature  end-heaters  were  sub¬ 
jected  to  a  constant  voltage  input  to  the  heaters  at  the  high  tem¬ 
perature  end  and  the  sink  allowed  to  achieve  thermal  equilibrium. 
The  heat  flowing  down  the  rod  vas  continually  removed  from  the 
low  temperature  end  of  the  calorimeter  by  iaserslng  that  end  in 
a  five  gallon  constant  te^e  rature  bath. 

F.  8Mjplea 

The  selection  of  the  test  materials  for  this  study  was  made 
on  the  basis  of  their  thermoelectric  properties.  Lead  vas  chosen 
for  the  Investigation  since  it  is  generally  accepted  as  a 
standard  thermoelectric  reference  material,  based  on  its  very  low 
Thomson  coefficient,  i.e,  almost  no  variation  of  thermoelectric 
power  vith  temperature  (cf.  Equation  l).  Thus,  any  thermoelectric 
data  obtained  between  lead  and  another  material  is  often  considered 
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to  to*  an  absolute  value  of  the  thermoelectric  properties  of  the 
second  material.  Vidcel,  on  the  other  hand,  vas  selected  on  the 
basis  of  its  vary  high  Thomson  coefficient  in  the  temperature  range 
to  be  studied.  The  choice  of  these  two  metals  should,  therefore, 
demonstrate  the  extreme  interactions  taking  place  in  say  measure* 
sent  of  the  Thomson  coefficient*,  the  difference  in  their  Thomson 
coefficient  values  is  approximately  two  orders  of  magnitude  * 

Pb  %  0.30  yv/deg.  and  Hi  %  35  yv/ deg. 

The  lead  used  in  this  investigation  vas  obtained  from  the 
American  Smelting  and  Research  Corporation  of  Rev  Jersey.  The 
major  constituents  in  the  lead  bars  vere  as  follows : 


Lead 

99.999  + 

Magnesium 

<  0.5  phi 

8111  con 

<  0.5  PPM 

Iron 

<  1.0  PPM 

Copper 

<  1.0  PPM 

Oxygen 

<  5.0  PPM 

Several  different  approaches  vere  used  in  an  attempt  to  ob¬ 
tain  a  0.020"  coaxial  hole  through  the  length  of  the  sample.  The 
most  satis  factory  results  vere  obtained  by  first  melting  the  lead 
in  a  pyrex  tube  under  an  argon  atmosphere  and  then  casting  the 
liquid  metal  into  a  5  mm  ID  pyrex  tube,  vhich  vas  held  at  325°C 
by  a  ni chrome  heating  coil  vrapped  around  the  glass.  Two  5  mm 
iron  slugs ,  with  25  mill  center  holes ,  vere  sealed  act  opposite  ends 

of  the  glass  to  hold  the  20  mill  iron  wire  in  the  center  of  the 
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tube.  A  three  pound  weight  vu  suspended  from  the  lower  end  of  the 
wire  to  hold  the  wire  straight .  To  prevent  leakage  of  the  liquid 
lead  out  of  the  wire  exit  hole,  wet  asbestos  paper  was  wrapped 
around  the  base  of  the  tube .  Once  the  lead  had  been  east  the  variac  . 
vas  gradually  reduced  to  zero.  The  glass  was  then  thermally  cracked 
and  finally  removed  from  the  specimen  by  careful  grinding.  Since 
iron  end  lead  are  completely  immiscible ,  the  iron  wire  could  usually 
be  removed  by  simply  pulling  the  wire  through  the  specimen. 

The  lead  sample  vas  then  mounted  in  the  calorimeter  by  placing 
one  end  of  the  sample  in  one  of  the  lead  heat  sinks  and  then  coating 
the  sample  where  it  met  the  sink  with  the  vacuum  varnish.  The 
sample  and  the  sink  were  then  carefully  guided  into  the  calorimeter. 
Once  the  sample  vas  in  place  in  the  opposite  end  of  the  calorimeter, 
the  sample-heat  sink  Junction  vas  given  a  light  coat  of  varnish  and 
the  second  lead  sink  slipped  into  place  over  the  sample  into  the 
copper  heating  sink.  At  this  time  the  thermocouple  vas  inserted 
and  the  leads  passed  through  the  electrical  connectors  as  described 
earlier. 

The  nickel  vas  obtained  from  the  United  Minerals  and  Chemical 
Corporation  of  Hev  York  City.  The  major  constituents  in  the  nickel 
rods  vers  as  follows i 
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Iron 


<  k  PPM 


8111 con  <  3  PPM 

Calcium  <  1  PPM 

Copper  <  1  PPM 

Magnesium  <  1  PPM 

Silver  <  1  PPM 

Since  it  vas  almost  impossible  to  cast  or  drill  a  0.5  to  1.0 
tm  hole  over  a  length  of  150  on  in  the  nickel  specimen,  the  sample 
had  to  he  prepared  from  tvo  5  am  diameter  rods.  Two  and  one-tenth 
centimeters  were  milled  longitudinally  from  one  of  the  nickel  rods 
and  then  an  .8  am  slot  vas  milled  the  entire  length  of  the  sample. 
Two  and  nine-tenths  centimeters  were  then  milled  from  the  second 
savle*  the  remaining  section  being  placed  on  top  of  the  first 
milled  section  and  the  two  section  then  Hell -Arc  welded  the  length 
of  the  rod. 

The  same  technique  vas  used  to  mount  the  nickel  specimen  in 
the  calorimeter  as  used  with  the  lead  specimen. 

G.  Experiment-*]  7*>lii)'in<,iTu*" 

There  were  basically  three  operational  techniques  employed 
in  obtaining  the  value  of  the  Thomson  coefficient.  Of  these 
methods,  tvo  were  essentially  new  techniques  developed  in  this 
investigation  from  the  fundamental  equations  governing  the  longi¬ 
tudinal  flow  of  heat  in  a  rod,  i.e. 

dx 

— Uo— 


(30) 


The  third  technique  vu  that  euggeeted  by  Nystrom,  in  vhlch  the 
Those  on  coefficient  vaa  calculated  froa  the  chaaga  in  current  nec- 
eaaary  to  compensate  for  the  change  in  heat  cfetalned  with  re  rare  al 
of  current  flow,  i.e. 

j2°  +  jT  S  "  (J  +  63)20  ’  (J  +  ^  S  (33) 

Method  I  ,  the  firat  technique  to  be  examined  vaa  baaed  on  non¬ 
equilibrium  considerations  of  Equation  (30).  Uhen  current  ia  first 
applied  to  a  simple  heater  coil,  almost  all  of  the  heat  ia  first 
absorbed  by  the  aaaple,  and  little,  if  any,  is  lost  to  the  surround¬ 
ings.  Mathematically,  this  is  based  on  the  equation  for  the  driving 
force  for  heat  losses: 

hA(Ts.  -  W  w 

Where 

h  ■  heat  transfer  coefficient 

*  temperature  of  the  surroundings 

Tga  ■  temperature  of  the  coil 

In  the  firat  fev  seconds  there  is  very  little  difference  between 
and  T&a,  and,  therefore,  thia  term  may  be  assumed  to  be  sero. 

Applying  thia  reasoning  to  the  basic  aquation  governing  heat 
flow  in  a  saaple  such  as  used  in  Thomson  heat  investigations,  the 
last  term  in  Equation  (30)  be  corns  s  sero,  and  the  equation  can  be 
rewritten: 


The  initial  runs  vere  mad*  on  the  pur*  lead  samples  (  in  vhlch 
the  current  vu  applied  in  one  direction  only  until  a  definite  elope, 
(dT/dt)^,  vac  established.  As  soon  as  this  slope  could  be  measured, 
the  current  was  then  reversed  and  a  second  slop*  established, 
(dT/dt)g,  due  to  the  different  rate  of  heat  input.  Mathematically, 
the  Thomson  coefficient  could  be  calculated  from  the  following 
considerations : 

For  current  in  one  direction: 

Bc(§>i  ■  &  ♦  Jt  S  (te) 

for  current  flow  in  the  opposite  direction: 

“<§>2  •  A  -  "  §  <«> 

subtracting  Equation  (43)  from  Equation  (42) 

acA  ft  ■  2Jt  fj  (44) 

Another  analysis,  based  on  the  same  non-equilibrium  equations, 
is  to  simply  measure  the  initial  slope,  (dT/dt),  for  either  direction 
of  current  flow  and  calculate  the  Thomson  coefficient  from  this 
measurement,  i.e. 

mc(f|)  -  A  +  Jx  f|  (45) 

T  *  7m  <“#  5  ‘ j2|,)  (W) 

dx 

or,  for  the  opposite  direction  of  flow: 

T  *  Tat lj2(>  ',c  £>  m 

dx 
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The  other  new  technique  used  to  calculate  t ,  Method  II,  Is 
hssed  on  equilibrium  considerations  of  the  basic  equation.  At  equi¬ 
librium  the  staple  is  not  changing  temperature  vith  time  and 
Equation  (30)  becomes: 

-  k  d^/dx2  -  J2p  *  Ji  ~ 

This  equation  has  been  solved  in  Appendix  B,  for  a  general 
position  (x)  in  the  sasple  and  then  specifically  for  x/2.  To  ob¬ 
tain  t  for  a  general  position  in  the  ss^l*  it  la  necessoxy  to  use 
a  computer.  The  known  parameters,  such  as  resistivity,  thermal 
conductivity,  end  temperatures ,  and  current,  were  substituted  into 
the  equation  for  various  values  of  t  .  The  counter  is  programed 

such  that  a  plot  of  both  directions  of  current  flow  was  made  of 
»• 

position  versus  temperature,  its  derivative  with  respect  to  po¬ 
sition,  and  the  difference  in  temperature  between  the  two  different 
directions  of  current  flow.  The  print  out  from  the  computer  for 
1  of  50  values  of  t  was  as  follows: 


x/L 

T_ 

6 

■  30  x  10 

6p 

■*  j. 

8' 

S 

46 

.1 

.1009 

1.008 

.1002 

1.001 

.0008 

.2 

.2017 

1.006 

.2003 

1.001 

.0014 

1.0 

1.0000 

•  9896 

.9999 

.9980 

.0000 
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where 


e 


8' 


dfi/dx 


for  current  in  one  direction 


for  current  in  opposite  direction 


By  substituting  actual  data  in  the  equation  and  comparing  the 
actual  experimental  values  with  the  set  of  values  obtained  from 
the  computer,  a  value  for  t  could  be  obtained. 

A  much  simpler  method  of  determining  t  is  to  measure  the 
change  in  temperature  at  the  center  of  the  saaple  and  simply  use 
the  equation  developed  in  the  latter  part  of  Appendix  B.  This  is 
by,  far  the  easiest  and  most  direct  method  of  measuring  the  Thomson 
coefficient.  It  has  been  developed  from  first  principles  and  com¬ 
pletely  does  away  with  any  need  of  changing  the  currant,  which  was 
found  to  be  rather  Insensitive  in  this  investigation. 


V.  RESULTS  AMD  DISCUSSION 


Method  I,  the  non -equilibrium,  change  in  slope  technique,  vas  used 
first  to  measure  the  Thomson  coefficient  of  lead.  A  multitude  of  runs 
were  made  with  these  samples ;  however,  stray  fields  were  of  such  magni¬ 
tude  ,  relative  to  the  small  changes  associated  with  the  Thomson  heat  of 
lead,  that  the  results  were,  in  general,  quite  erratic.  The  only  re¬ 
producible  results  were  obtained  late  at  night  when  the  fields  around 
the  measuring  equipment  were  at  a  minimum.  Two  of  the  more  stable  runs 
on  the  lead  sasples  are  shown  in  Tables  1A  and  IB. 

TABU!  1A 


Values  of  the  Thomson  coefficient  of  Lead  Using  Method  I 


Base  Temperature  ■  Ul.0°C 

dET/dx  “  2.2°C/cm 

Amps 

dpv/dt 

Thomson  coefficient 
jit  volts /degree 

3.0 

0.0008 

-  .412 

3.96 

0.0009 

-  .339 

5.10 

0.0013 

-  .380 

5.98 

0.0016 

-  .398 

6.85 

0.0013 

-  .283 

Average  t  ■  0.360  *  O.Obl 
x  10“^  volts /deg. 

TABLE  IB 


Values  of  the  Thomson  coefficient  of  Lead  Using  Method  I 


Amps 

dyv/dt 

Thomson  coefficient 
p  volts /degree 

3.92 

0.0005 

-  0.189 

5.16 

0.0017 

-  0.498 

5.98 

0.0009 

-  .235 

6.67 

0.0013 

-  .283 

Average  0.304  *  0.099  x  10  “6 
volts/degree 

One  of  the  calculations  is  as  follows: 


for  lead 


me  dT/dt  «  2Jt  dT/dx 

m  »  11  ©n/cc  (44) 
c  ■  .03  cal/©a-°C  (44) 


(n  ®s/cc)(.03  -£sji)( — i-sa£t — >/ .  2  (3.iP,«gg)(2>1 

J  @b-®Cm.239  cal/secM42.1mr;v  dt'  *  '  2  n*,A  os' 


t  *  .412  x  10"°  pv/deg 


The  average  results  obtained  with  this  technique  is  fairly  close 
to  the  value  of  .30  yv/ deg,  given  by  Christian  et  al  (45);  however,  the 
scatter  is  fairly  high.  To  test  the  validity  of  this  experimental 
approach,  two  alterations  were  made  in  the  approach.  The  first  change 
was  made  in  order  to  make  the  measurements  further  from  thermal  equi¬ 
librium,  where  the  heat  losses  were  more  likely  to  be  sero.  This  was 

accomplished  by  changing  to  much  faster  recording  rates  in  the  Brown 
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Electronic  Recorder  and  using  a  single  slope  to  calculate  the  Thomson 
coefficient,  i.e. 


dT  _2  A  T  <1T 
me  —  -  1  p  +  Tt  — 


The  second  change  was  to  simply  use  nickel  samples,  where  the  Thomson 
coefficient  is  so  high  that  the  stray  fields  would  only  be  a  small  per¬ 
centage  of  the  Thomson  heat. 

The  experimental  results  obtained  using  this  technique  were  found 

to  have  much  less  scatter  when  the  electron  flow  was  in  the  direction 

of  decreasing  temperature .  The  slopes  of  the  curves  for  this  direction 

of  flow  were  quite  linear  in  the  first  few  seconds  and  gradually  changed 

slope  as  the  heat  losses  became  more  pronounced.  With  electron  flow  in 

the  opposite  direction,  i.e.  in  the  direction  of  increasing  tespera- 
•> 

ture,  it  was  difficult  to  obtain  reproducible  results.  During  the  first 
few  seconds  of  such  a  run,  a  definite  slope  would’  first  be  established, 
but  would  suddenly  give  way  to  a  much  steeper  slope  as  shown  in  Figure 
6.  This  behavior  is  believed  to  be  due  to  sosk  fundamental  process  tak¬ 
ing  place  in  the  saople,  which  is  not  encountered  in  the  more  conaon 
equilibrium  considerations.  In  making  calculations  of  the  Thomson  co¬ 
efficient  with  opposing  thermal  and  electron  flow,  the  earliest  slope 
was  used  to  make  the  calculations.  This  is  admittedly  questionable, 
but  the  results  do  agree  better  with  those  obtained  with  the  opposite 
flow,  although  they  are  not  nearly  as  reproducible. 

The  average  of  four  runs,  ignoring  small  differences  in  experi¬ 
mental  conditions,  is  shown  in  Table  2. 
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Temperature,  °C  Temperature,  *C 


Case  I  *  (Electron  Flow  Down  Thermal  Gradient) 


.1 


05 


/ 

- 1 - 4— 

5  10 


Time,  sec. 


Time,  sec. 


FIOURS  6  -  Sketch  of  curves  obtained  with  two  directions  of  current 
flow  in  a  nickel  sample 
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TABLE  2 


Values  of  the  Thomson  coefficient  of  Nickel,  Using  Method  I 


TBase  "  hl •2°C  "  2*5°C 

Amps 

A* 

dT/dt,°C/sec  -Ttuv/deg 

B* 

dT/dt ,°C/sec  -t ,uv/deg 

2.98 

3.96 

5.08 

5.96  ' 
6.78 

.0094  24.74 

.0197  30.93 

.0279  39.74 

.0464  30.79 

.0710  35.52 

.0318  42.8 

.0381  13.81 

.0676  36.60 

.0822  20.05 

.1340  45.20 

Average  —32.34  *  4.23  vv/deg 

Average  *-31.69  *  10 .81 
ur/deg 

Where  A*  ■  values  for  electron  flow  in  the  direction  of  decreasing 
„  temperature 

B*  ■  values  for  electron  flow  in  the  direction  of  increasing 
temperature 

A  sample  calculation  is  as  follows: 

I  *  3.96  amps 


where  for  nickel 

m  ■  8,75  ffn/cc  (46) 
c  ■  .11  cal/gm-°C  (46) 


(6.75H.LL)  £U 


(7.7)  *  10"*  -  (2.5)t 


t  ■  30.93  w/dag 


-49- 


Several  runs  were  made  using  the  standard  technique  of  changing  the 
current  to  condensate  for  the  difference  in  the  heat  produced  on  switch¬ 
ing  the  direction  of  current  to  calculate  the  Thomson  coefficient,  i.e. 
Method  III.  This  technique  was  used  to  see  how  well  the  experimental 
apparatus  was  performing  with  respect  to  the  arrangement  used  by  other 
investigators.  The  exact  technique  employed  for  this  measurement  can 
best  be  understood  with  reference  to  Figure  7.  The  current  was  first 
applied  to  the  sample  such  that  the  Thomson  heat  was  adding  to  the  Joule 
heat  and  the  thermocouple  voltage  noted  when  equilibrium  was  obtained, 
point  1  in  Figure  7*  The  current  was  then  reversed  and  the  equilibrium 
voltage  measured  with  this  direction  of  current  flow,  point  2.  The 
current,  with  the  direction  of  current  now  opposite  to  that  first  ap¬ 
plied,  was  then  arbitrarily  increased  and  the  associated  change  of  the 
thermocouple  was  measured  once  an  equilibrium  temperature  was  achieved, 
point  3.  The  current  was  further  adjusted,  points  U  and  5 ,  until  the 
equilibrium  plateau  was  very  close  to  that  first  achieved,  point  1. 

This  data  could  be  used  in  a  direct  calculation  of  the  Thomsen  coef¬ 
ficient,  from  Equation  (33);  however,  a  better  average  result  can  be 
obtained  for  a  given  set  of  boundary  conditions  by  continuing  the  sane 
procedure  to  higher  currents.  Starting  at  point  5,  the  current  can  be 
increased  to  some  new  value,  point  6,  noting  associated  changes  in  enf 
output,  and  then  reversing  the  current  flow,  point  7.  The  current  could 
then  be  decreased  arbitrarily,  point  8,  and  then  more  exactly,  points  9 
and  10,  until  the  emf  corresponding  to  point  6  is  obtained*  With  the 

experimental  arrangement  used  in  this  investigation  the  current  was 
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fflci tat 


i 

i 
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vailed  up  to  approximately  8  asps  and  then  bach  down  to  3  asps . 

In  order  to  cover  this  current  range  and  maintain  the  precision 
necessary  for  these  calculations,  it  vas  necessary  to  cancel  part  of 
the  signal  from  the  thermocouple  in  going  from  one  measurement  of  the 
Thomson  coefficient  to  another,  i.e.  from  points  3  and  6.  The  read¬ 
justment  might  possibly  have  introduced  a  slight  shift  in  the  data  be¬ 
tween  points  such  as  5  and  6.  It  should  be  emphasized  that  this  in  no 
way  affects  the  calculated  values  of  the  Thomson  coefficient,  since  no 
adjustments  were  necessary  in  the  series  of  points  1  through  5  and  6 
through  10.  Although  there  might  be  some  small  errors  by  presenting 
the  data  as  a  graph ,  it  was  felt  that  this  type  of  presentation  pro¬ 
vided  much  better  means  of  comparison  with  data  calculated  by  other 
means.  ~  Figures  8,  9,  and  10,  are  the  graphs  plotted  from  the  data  ob¬ 
tained  in  this  way.  The  resulting  values  of  the  Thomson  coefficient, 
calculated  from  the  equation: 


T 


pAJ 

ar/dx 


are  shown  on  the  respective  graphs .  The  dashed  lines  in  Figure  10  is 

a  plot  of  the  results  obtained  from  the  computer  technique  using  values 

of  t  indicated  on  the  graph,  where  the  boundary  conditions  are  the  same 

as  those  using  the  experimental  run.  The  results  obtained  from  the  two 

different  techniques  appear  to  agree  reasonably  well. 

Calculations  using  the  simplest  of  all  techniques ,  developed  in 

the  latter  part  of  Appendix  B,  are  shown  in  Table  3.  This  technique 

could  be  applied  to  all  the  data  obtained  in  Method  III,  since  it  is 
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TB«m  ■  *‘0-0“0  »  "l  *  26-9°C  1  **  •  51.2“C 


FIGURE  8  -  Data  obtained  on  nickel  samples  using  Method  III  - 
specific  boundary  conditions  (numbers  in  body  of 
graph  are  values  of  dl) 
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TBase  “  l‘2,15°c  *  T2  ■  55.1°C  ;  ?1  ■  30.0°C 


FIGURE  9  -  Data  obtained  on  nickel  samples  using  Method  III 
specific  boundary  conditions  (numbers  in  body  of 
graph  are  values  of  AI) 
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also  an  equilibrium  technique  using  data  obtained  at  x/2.  The  only  dif¬ 
ference  being,  that  this  latter  technique  is  much  easier,  in  that  the 
current  need  not  be  changed.  The  only  data  necesaaxy  ia  the  change  in 
equilibrium  temperature,  which  is  the  first  data  obtained  in  Method  III. 
The  Thomson  coefficient  can  simply  be  calculated  from  the  equation: 


T 


kk.LT 

W1 


(M) 


TABLE  3 


Values  of  the  Thomson  coefficient  Using  Method  II,  Special  Case 


V 

■  27.8°C  T2  ■  52.6°C  TCm<j  -  40.1°C 

k  -  .594  watts/cm  (44) 

L  *  9.3  cm 

Amps 

AT 

-r,  yv/  degree 

2.96 

CO 

j* 

o 

e 

33.3 

4.02 

.048 

24.6 

4.88 

.067 

28.2 

5.98 

.107 

36.8 

6.92 

.098 

29.0 

Average  *-32.4 

*  4.0  x  yv/deg 

For  the  sake  of  comparison ,  this  equation  has  been  used  to  calcu¬ 
late  the  values  of  t  from  the  data  used  in  Figure  10.  The  values  of 
the  Thomson  coefficient  obtained  in  this  way  are  given  in  Table  4.  Any 
difference  between  the  2  techniques  is  probably  due  to  small  radial 

heat  losses  which  were  not  taken  into  account  in  the  equation  developed 
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in  Appendix  B  and  the  fact  that  different  coefficients  vers  used  in  the  . 
two  techniques  to  calculate  the  Thomson  coefficient. 

TABLE  4 

Thomson  coefficient  Calculations  from  Previous  Data,  Using  Method  II  - 

Special  Case 


Amps 

dyv 

-r.yv/deg 

3. 

1.81 

29-9 

3.5 

2.17 

31.0 

4.0 

2.50 

31.2 

4.5 

2.85 

31.6 

3.22 

32.2 

3.57 

32.5 

3.95 

32.9 

4.25 

32.7 

Average  *-31.8 
yv/  deg 

Values  of  the  Thomson  coefficient  of  nickel  and  lead  in  the  tem¬ 
perature  range  studied  in  this  investigation  are  quite  limited.  In 
order  to  make  a  direct  comparison  of  these  data  with  the  results  ob¬ 
tained  from  other  studies ,  it  vaa  necessary  to  make  extrapolations  of 
both  Thomson  coefficient  data  and  thermoelectric  power  data,  which  was 
then  differentiated  to  obtain  values  for  the  Thomson  coefficient. 

Using  the  results  obtained  by  Bo  re  li  us  (47),  the  Thomson  coef¬ 
ficient  of  lead  in  the  tesgmrature  range  of  40.0°C,  would  be  -  0.57 
yv/ deg,  which  is  in  rough  agreement  with  the  average  value  of  -  0.33 

*  .09  yv/de  g  found  in  this  study.  Bore  li  us  indicated  that  his  values 
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could  have  been  in  error,  as  his  equipment  had  not  been  properly  stan¬ 
dardised.  Extrapolating  the  thermoelectric  power  results  obtained  by 
Pearson  et  al  (40),  to  the  teiqperature  range  around  4o.O°C,  and  then  tak¬ 
ing  the  derivative  of  this  curve,  values  of  approximately  -  0.30  yv/ deg 
are  obtained  for  lead,  which  agrees  quite  closely  with  the  values  obtained 
for  lead  using  the  change  in  slope  technique.  The  values  calculated  with 
the  change  in  slope  technique,  however,  are  quite  scattered.  Since 
Pearson  has  not  reported  any  limits  in  his  work,  it  is  difficult  to  make 
a  comparison  between  the  two  means  of  measurement. 

Hi eke 1  has  been  studied  by  Borelius  (U7) ,  Budnitskii  (13)  ,  and 
Grew  (48).  Borelius*  results  indicate  a  value  of  -  18.0  vr/ deg  at  4l.0°C, 
although  he  freely  admits  that  this  value  is  probably  not  too  valid, 
since  hit  samples  were  relatively  impure  and  the  calculations  were  made 
with  a  questionable  standard.  Budnitskii's  results  ihdicate  at  value  of 
approximately  -  32.0  yv/deg  for  nickel  at  Ul.0°C,  which  agree  much  better 
with  the  value  of  approximately  -  30.0  *  $.0  yv/deg  obtained  in  this  in¬ 
vestigation.  He  has  not,  however,  given  any  indication  of  the  amount  of 
scatter  which  may  have  been  present.  Extrapolating  Grew's  measurements 
of  thermoelectric  power  to  Ul.0°C  and  taking  the  derivative  of  this  curve, 
yields  a  value  of  approximately  -  28.0  yv/deg  for  nickel,  which  again 
agrees  very  well  with  the  data  obtained  in  this  study.  Once  again,  no 
indication  is  given  in  regard  to  the  possible  amount  of  error. 

As  indicated  above ,  the  values  of  the  Thomson  coefficient  of  both 
lead  and  nickel  obtained  in  this  investigation,  were  apparently  the 
first  to  be  calculated  from  direct  measurements  in  the  temperature  rsnge 
of  40.0°C. 
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There  seem*  to  he  no  apparent  explanation  for  the  change  in  elope 
encountered  with  opposing  thermal  and  electrical  flows ,  in  the  non -equi¬ 
librium  technique.  If  this  phenomena  were  due  to  soma  type  of  phonon- 
electron  interaction,  one  would  expect  the  effects  to  be  completed  in  a 
matter  of  micro-seconds.  As  this  anomaly  was  generally  encountered  a 
second  or  two  after  the  initiation  of  current  flow,  such  an  explanation 
seems  quite  inadequate.  Magnetic  effects  associated  with  the  initiation 
of  the  D-C  current  were  also  ruled  out,  since  these  fields  would  be 
stationary  in  a  very  short  time.  To  insure  that  effects  such  as  these 
were  not  delayed  between  their  generation  and  recording,  a  rapidly  chang¬ 
ing  voltage  from  the  Lin  deck  microvolts  source  was  impressed  on  the  mea¬ 
suring  network ,  which  indicated  that  the  circuit  had  a  very  fast  response 
time.  This  phenomenon  does  seem  to  be  attributable  to  some  interaction 
taking  place  in  the  sample,  which  only  occurs  with  opposing  thermal  and 
electrical  flows. 
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VI.  SUMMARY  AND  RECOMMENDATIONS 


In  general ,  this  investigation  presented  the  fundamental  concepts 
or  the  Thomson  heat  and  applied  these  concepts  to  measure  values  of  the 
Ttmomson  coefficient.  More  specifically,  the  investigation  can  be 
lesparated  to  include  the  following: 

1.  Development  of  the  non-equilibrium  techniques  for  the  direct 
measurement  of  the  Thomson  coefficient. 

a)  This  included  the  change  in  slope  technique,  and 

b)  Single  slope  technique. 

2.  Development  of  a  nev  approach  to  the  equilibrium  techniques 
for  determination  of  the  Thomson  coefficient. 

a)  This  included  a  method  for  a  general  position  in  the 

w 

sample,  and 

b)  The  special  case  of  measurements  made  at  x/2. 

3.  The  design  and  utilization  of  an  experimental  apparatus  to 
examine  the  various  approaches  for  a  direct  determination  of 
the  Thomson  coefficient.  The  samples  in  this  investigation 

s 

vere  designed  to  avoid  erroneous  temperature  measurements. 

In  either  welding  the  thermocouple  to  the  side  of  the  specimen 
or  drilling  surface  holes  In  the  specimen,  the  actual  tem¬ 
perature  distribution  and  electrical  flow  is  disturbed  in  the 
area  at  which  the  measurement  is  being  made.  Although  the 
affect  may  be  rather  small,  it  could  amount  to  an  appreciable 


error  in  relation  to  the  small  changss  associated 
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with  the  Thomson  heat.  With  the  placement  of  the  thermocouple 
down  the  center  of  the  sample  t  as  in  this  investigation ,  these 
effects  would  be  minimized;  this,  however,  may  introduce  small 
losses  of  heat  down  the  center  of  the  sasple,  which  undoubtedly 
have  some  small  affect  on  the  experimental  values.  Theoreti¬ 
cally,  the  beat  method  for  avoiding  theBe  difficulties  is  the 
approach  used  by  Nystrom,  i.e.  an  optical  thermal  measurement. 

In  practical  application  this  approach  has  proved  to  be  no 
more  reliable  than  the  others ,  since  it  is  usually  limited  by 
endssivity  of  the  radiant  energy  from  the  sample  and  absorption 
in  the  optical  path. 

Application  of  these  new  techniques ,  as  well  as  standard  tech¬ 
niques  ,  to  calculate  the  Thomson  coefficient  of  lead  and  nickel. 

a)  Method  I  -  Hon -equilibrium  Technique  -  This  technique 
proved  to  be  fairly  reliable,  with  a  minimum  of  scatter, 
for  electron  flow  in  the  direction  of  decreasing  tempera¬ 
ture.  The  results,  however,  were  not  reproducible  for  the 
current  flow  in  the  opposite  direction. 

b)  Method  II  -  Hew  Equilibrium  Techniques  -  The  special  case 
of  this  technique,  i.e.  measurements  made  at  a  point  half 
way  between  the  two  end  sinks,  is  by  far  the  simplest 
approach  for  calculating  the  Thomson  coefficient.  This 
technique  does  require  very  small  radial  heat  losses. 

c)  Method  III  -  Standard  Technique  -  The  adjustment  of  the 

current  to  exactly  compensate  for  the  loss  or  gain  of 
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Thomson  heat  was  found  to  he  quite  difficult  to  control 
with  the  experimental  arrangement  used  in  thie  investi¬ 
gation.  The  results,  however ,  were  fairly  reproducible 
and  reliable,  as  long  as  the  assumption  is  made  that 
there  is  no  change  in  the  heat  loss  mechanism  with  change 
in  current  flow. 

The  measurement  of  the  small  changes  occurring  in  this  investi¬ 
gation  were  hampered  by  stray  electromagnetic  fields  in  the  laboratory. 
This  was  validated  by  purposely  introducing  electromagnetic  source  or 
magnetic  fields  during  an  experimental  run.  The  stray  fields  were  cir¬ 
cumvented  in  part  by  shielding  the  lead  vires  as  much  as  possible,  by 
minimizing  loops,  and  making  runs  late  at  night  when  the  fields  were  at 
a  minimum.  It  is  felt,  however,  that  the  accuracy  and  reproducibility 
of  the  data  was  limited  by  these  fields  .  For  this  reason  it  is 
recommended  that  future  studies  of  this  nature  be  conducted  in  a  well 
grounded  cage. 

The  anomaly  encountered  with  opposing  electri  cal  and  thermal 
currents  presents  an  interesting  problem,  the  solution  of  which  was  be¬ 
yond  the  scope  of  this  investigation.  An  investigation  of  the  mech¬ 
anisms  giving  rise  to  such  an  effect  could  possibly  lead  to  a  better 
understanding  of  electron-phonon  interaction  and  should  be  pursued  on 
this  basis. 
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APPENDIX  A 


A  simple  heater  coil  will  be  considered  to  show  hew  various  equi¬ 
librium  temperatures  can  be  achieved  by  simply  changing  the  rate  of 
heat  input.  The  general  equation  governing  the  flow  of  heat  into  the 
coil  is  as  follows: 

Qj  -  Qa  ♦  Qj,  (1A) 


or  more  explicitly 


I‘pv  -  cm  ♦  hAtT^  -  TSu) 


(2A) 


where: 


I  pV  ■  rate  of  electrical  heat  input 


cm 


dT 

dt 


“'VW 


rate  of  heat  accumulated  in  the  coil: 
c  -  specific  heat  and  a  «  density 

rate  of  heat  losses  to  the  surrounding: 

h  ■  heat  transfer  coefficient 

A  ■  area 


Tga  ■  sample  temperature 


Su 


temperature  of  surromdings 


This  equation  has  a  solution  of  the  form: 

t2  v  T2flV  -te 
<r  „  1  PV  1  PVe _ +  T 

Sa  hA  hA  Su 


A  plot  of  the  various  heats  versus  tine  is  shown  in  Figure  11. 

According  to  Equation  (IB),  if  the  rate  of  heat  input  is  now  in¬ 
creased  the  previous  curves  will  be  shifted  up  as  shown  in  Figure  12, 
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coning  to  a  different  equilibrium  temperature ,  assuming  the  heat  transfer 
coefficient  does  not  change  with  change  of  current  flow. 

e 

Inclusion  of  the  Thomson  heat  in  the  term  in  the  previous  analysis 
will  not  have  an  appreciable  affect  on  the  shape  of  the  curves . 


APPENDIX  B 


The  •inpHfied  equilibrium  equation  governing  heat  flow  out  of  a 
rod,  vhere  radial  heat  loaaea  are  considered  to  he  negligible,  is: 

-  kd^/dx2  M  J2/o  +  tJ  (IB) 

°*1 

The  solution  to  this  equation,  in  dimensionless  analysis  to 
facilitate  use  in  the  computer,  is  as  follows: 
assuming  the  boundary  conditions , 


T(0)  ■  T  i  T(L)  ■  T. 

o  u 


and  defining  a  dimensionless  temperature  as, 

8  • 


»<x>  •  To 


TL-To 


where  ■  the  temperature  at  some  point  (x)  along  the  ssmple  between 
0  and  L,  such  that, 

8(0)  •  0  i  8(L)  -  1  ; 


and  letting 


then 


Letting 


and 


x  ■  x^L, 

d^e.  Jtl  de  .  j^l2 
ta2  k  ax  k«(T^o) 

a  ■  IxL/k 

b  -  I2L2/ko(T.  -  T  ) 
Xj  O 


(2B) 
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Equation  (2B)  becomes: 


a 


de. 

dx 


+  b  ■  0 


Let  0  ■  0^  + 


6p  where  is 


the  solution  to  the  homogeneous 


equation: 


0 


let 


d6h  dn 

p  then*  -  «P 


/ 


/  a  dx 


P 


*  C  e 

o 


-ax 


CQ  *  integration  constant 


!V.c  .  “ 

dx  o 


—  e  +  C 


6p  is  a  particular  solution  given  by 

Vx) " » x 


e---2-  ♦  c, 

*  ft  JL 


c  c 

e(o)  ■  0  ■  -  —  e”a^  +  c.  .*•  c..  -  ~ 

ft  <L  X  ft 
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c  .  c  c 


e(i)  -  1  -  (!-«•*)  -  £ 


a+b 


(l-e'a) 


a  (l-e"a)  a 


(3B) 


For  current  flew  in  the  reverse  direction: 


d2e„  dflp 

- 1  -  a  -~+  b  •  0 

dx2  ** 


therefore,  to  solve  the  reversal  current  flow  problem,  substitute 
a  ■  -a,  and 

eRU)  "  '  +  ~  *  (to) 

R  a  (l-ea)  a 

Equations  (36)  and  (46)  are  the  general  solutions  to  the  equation 
governing  longitudinal  heat  flow  in  the  rod.  Calculations,  using  the 
previous  equations,  of  position  versus  tenperature  over  a  linear  base 
are  shown  in  Figure  13.  The  same  equations  predict  the  slope  to  change 
in  the  saople  as  shown  in  Figure  l4. 

It  is  interesting  to  consider  the  change  in  equilibrium  tempera¬ 
ture  at  the  middle  of  the  rod  when  the  current  is  reversed. 

0  -  6-  ■  A8  ■  +  t,aZ&l  .  . )  _  (5B) 

R  »  (1-.-*)  *  <!-.»)  * 
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It  should  he  noted  that  this  result  is  independent  of  the  para¬ 
meter  h,  which  includes  the  Joule  heat  end  the  electrical  conductivity 
of  the  material. 

A  series  expansion  of  the  previous  exponentials  gives  a  very 
simple  expression  fpr  Equation  (5B). 


Expressing 


n  2n  2n+l 

*  Jo  Jr  ■  Jo  7&TT*  Jo  TfciJT 


and 


2n  2n+l 

«"a  "  nlo  ifen  *  n*0  TfelTT 


^m/O  a/2 

and  similarly  for  e  and  e  .  Substituting  in  Equation  (5B)  and 
siqpli  lying 


2n+l 


48 


Jl  ^  ^  (  2n+l)! 

I  s 

nil  2n7 


Using  the  first  tvo  terns  of  each  series  yields: 


48  ■ 


2  k 

t*  t 


or 


46  -  0.99999  x  £ 


46  ■ 


T<»> ~  V  . 

T.  -  T  F 

L  O 

1T  -  (IL  -  To)  f 

«  *  «L  ’  V  TT 


T 


Uk4T 

\  -  T0)Jl 
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APPENDIX  C 


List  of  Symbols 


V 

ir 

T 

T 


Seebeck  potential,  volune 
Peltier  heat 
Thomson  coefficient 
temperature 


maximum  temperature 
temperature  of  sample 

temperature  of  surroundings 

temperature  at  end  of  sample 

temperature  at  opposite  end  of  sample 

temperature  at  intermediate  point  in  sample 


e 

t 

s 


dimensionless  temperature 
time 


T 

? 


jx)  ~  T(0> 
<L)  *  T(0) 


thermoelectric  power,  entropy 
entropy  density 


6 

« 

Q 

cg 


entropy  flow 

transport  entropy  per  particle 
transport  heat  per  partirle 
lattice  specific  beat 


cel 

a 

8 


electronic  specific  heat 

probability  of  phonon -electron  collision 

number  of  conduction  electrons  per  atom 
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A 

A(B) 

t 

E 


« 


J 

e 


M  L 
eq  e* 


U 

y 

I 

J 


P 

R 

a 

k 


▼ 

x.y.z 

f 

h 


charge  on  electron,  logarithmic  base 
cross  section 

energy  dependent  scattering  cross  section 
mean  free  path 
electron  energy 

energy  of  electron  at  Fend  surface  at  0°K 

electric  and  heat  current  flow  vectors 

forces  giving  rise  to  current  flows 

phenomenological  coefficients 

particle  current  density 

entropy  flow  density 

interaction  coefficients 

chemical  potential 
electrochemical  potential 
electrical  current 
electrical  current  density 
resistivity 
resistance 

electrical  conductivity 

thermal  conductivity,  Boltzmann  constant 

velocity  of  electrons  in  specific  direction 

distribution  function 

Plank's  constant,  heat  transfer  coefficient 
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*  electron  density 

*  Fend  energy  per  electros  jsr  electron  charge 

■  rate  of  heat  input 

*  rate  of  heat  accumulated  tey  a  aaple 

*  rate  of  heat  leaving  a  unit  voluan 

*  electrical  potential 

■  specific  heat 

■  length  of  saxaple 
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